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ABSTRACT 


In the presence of 3-D turbulence, peak horizontal velocity estimates from an 
idealized Doppler profiler are found to be positively biased due to an incomplete 
specification of the vertical velocity field. The magnitude of the bias was estimated by 
assuming that the vertical and horizontal velocities can be separated into average and 
perturbation values and that the vertical and horizontal velocity perturbations are 
normally distributed. Under these assumptions, properties of the Type-I Extreme Value 
Distribution for maxima, known as the Gumbel distribution, can be used to obtain an 
analytical solution of the bias. The bias depends on geometric properties of the profiler 
configuration, the variance in the horizontal velocity, and the unresolved variance in the 
vertical velocity. When these variances are normalized by the average horizontal 
velocity, the bias can be mapped as a simple function of the normalized variances. 
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1. Introduction 


Doppler profilers are used to obtain estimates of horizontal ( U ) and vertical 
velocity ( W ) within the atmosphere (Van Zandt 2000) and under water (Woodward and 
Appell 1986) using acoustic, radar and optical remote sensing techniques. 
Electromagnetic or acoustic energy at a known frequency is transmitted into the medium 
of interest and the frequency of the back-scattered energy is measured by a directional 
receiver. Receiver characteristics such as size and shape define its beam, typically a 
narrow cone projecting away from the receiver throughout some depth of the fluid. The 
difference between the transmitted and received frequencies, referred to as the Doppler 
shift, is used to estimate the velocity component of back-scatterers along the beam axis. 
The back-scattering elements are assumed to be passive tracers of the fluid motion and 
the estimated velocity component along the beam axis is referred to as the radial velocity 
(R). Sophisticated transmitter/receiver configurations and signal generating/processing 
techniques have been developed over the past several decades to maximize the accuracy 
of R estimates and subsequent retrieved U and W . 

The average retrieved U (t/ Re , ) from Doppler profiling systems has been validated 
extensively with in situ observing systems in both the atmosphere (May et al. 1989 ; 
Crescenti 1997) and underwater (Chereskin et al. 1987; Gilboy et al. 2000). U Rel is 

generally well-estimated for averaging times of 15 minutes or greater. Average measures 
of turbulence can also be derived from Doppler profiler observations by 
statistical/dynamical methods (Kramar and Kouznetsov 2002; Lu and Lueck 1999). 
However, the precision of the retrieved instantaneous U ((t/ + u ) Ref ), where u implies 
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a perturbation from U with a time scale of a few seconds, is more problematic. As a 
result, the accuracy of the retrieved peak U ( {u® ) Re( ) from a collection of (u + w’)^, of 
size n, can be significantly affected. 

This paper describes an idealized Doppler profiler in an idealized fluid, where the 
true instantaneous R ((r +r') Tnie ) is composed of weighted sums of (u +u') True and 

(w + w ) 7>ae , with weights dependent on the beam configuration. In the case of 3-D 
turbulent flow the typical beam configuration does not provide adequate information for 
an accurate estimate of (u +u ) Re , , although ( U +u ) Re , can be shown to be unbiased 
under certain assumptions. For applications where U® from a collection of U + u of 
size n is of interest, errors in (u + w') Re( propagate into \U® ] Rel and, in general, introduce 

a positive bias in the average peak value {U® ) Re( . Section 2 describes an idealized 
profiler configuration along with retrieval algorithms for uniform and turbulent flows. 
Section 3 uses Extreme Value Theory to provide analytical solutions to (u® )^ e( and 

as a function of hirbulent properties of the fluid and the profiler configuration. 
Section 4 provides a summary and conclusions. A list of symbols is given in appendix A. 

2. An Idealized Doppler Profiler 

The following description of an idealized Doppler profiler is intended to represent, 
in the simplest terms, how {u +u') True and (w +w') Tnje combine to form (r + r) Tnie from 
the oblique and vertical beams of a typical system. The {r + r ) rree are then used to 
obtain {u +u ^ e( . For conditions where the (w +w ) Tnje 
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varies in space and/or time, 



(y+o'L is shown to be susceptible to an error that is directly proportional to the 
difference in ( W + w ) Tne over the oblique and vertical beams. A statistical modeling 
approach is used to show how this difference affects (u® ] Re; . The statistical model 
quantifies three important characteristics: 1) (U + w) Re , is unbiased under reasonable 

assumptions, 2) (t/® ) Rel is positively biased under the same assumptions, and 3) The 
magnitude of the positive bias is dependent on the unresolved temporal and/or spatial 
variations of ( W + w ) True , expressed in w 2 -w\. 

Consider an idealized Doppler profiler that measures the ( R +r ) Tne along each of 
two beams, bi and hi, as in Figure 1. The bi-beam is vertically oriented with respect to 
the local horizontal plane. The b 2 -beam is oriented at an angle 0 from the vertical in 
order to obtain information about U . While the typical Doppler profiling system has 
three or more beams to resolve the three orthogonal velocity components, the essence of 

the mathematical and statistical arguments supporting a positive bias in (t/®)^ e( can be 
readily developed for a two-beam system and are generally applicable for a multi-beam 
system. 

The profiler obtains doublets of (r +r ) Tne , Vj (Dj) and V 2 (D 2 ), where the 
symbols Dj and D 2 denote distances from point p in Figure 1, the location of the 
instrument. The condition Di = D2-Cos(0) is required to obtain radial velocities from 
both beams at the same height, H = Di, above the local horizontal plane. The following 

r 

notation was developed by assuming this condition was met and dropping the notation for 
height and distance. Note also that 0 is typically about 15° for atmospheric systems and 
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30° for underwater systems in order to obtain profile information over a useful depth and 
to avoid \>2 side-lobe contamination from ground targets with zero-Doppler shift. 
Because © is small, V 2 can be significantly affected by (w + w\ me . 

a. Equations for a Uniform Velocity Field with no Turbulence 

Under idealized conditions of temporal and spatial uniformity in the velocity field, 
the idealized profiler observations are described by the following equations (Peterson 
1988): 

V. = (W'L, (1) 

V 2 = (U) T „, ■ sin(0) + (w) T „, • cos(0) (2) 

Equation 2 can be directly solved for (U ) True as follows: 

TV, “ V 2 • sec(0) - (WV. ■ «>t(0> (3) 

There are three important points to note from Equations 1-3: 

1. Equation 1 shows that (W ) True is obtained from Vj, the radial velocity 
measured by the vertically oriented bi-beam. 

2. (W ) True appears in the second term for the (U ) True solution in Equation 3. The 
second term makes a correction for the effect of (W ) Tnje on V 2 . 

3. The correction term is amplified by the cotangent of the oblique beam angle, 
0. For 0 = 15°, cot (15°) = 3.73 and sec (15°) = 3.86. For 0 = 30 °, cot (30 °) 
= 1.73 and sec (30 °) = 2. 
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b. Equations for a Turbulent Velocity Field 

Consider a turbulent fluid where the (w) Tnie and (U ) True can be expressed in terms of 

average ( ) and perturbation ( ' ) values, (iV + w') True over the bi- and b 2 -beams will be 
denoted by Wi,= W t +wj and W 2 . = W 2 +w 2 , respectively, where W,=W 2 = 0. The 
(if + u ) True over the b 2 -beam will be denoted by U 2 +u 2 , where U 2 >0 will be assumed. 
In order to characterize a retrieval algorithm for these turbulent conditions, consider the 
following revised formulation of the profiler observations: 

Vi = wj , (4) 

V 2 = ip 2 + u 2 ) True • sin(0) + w 2 • cos(0). (5) 

The appropriate solution for the true horizontal velocity would be 

ip 2 +U 2 ) T rue = V 2 * Sec(0) - W 2 • COt(0). (6) 

However, in general, w 2 is not observed. It is approximated by w,, as in the 
following equation for the retrieved horizontal velocity: 

(U 2 + u 2 ) Re( = V 2 • sec(0) - w[ • cot(©) (7) 

(U 2 + u 2 ) Rer can be expressed as the true horizontal velocity, (ll 2 + u 2 ) Tnie , plus an 
error term, by combining Equations 6 and 7. 

(u 2 +u 2 ) Re , = (U 2 + u 2 ) True +(w 2 -w\)- cot(0) (8) 

The error term in Equation 8 is composed of the difference between the perturbation 
vertical velocities from the two beams, amplified by the cotangent (0) factor. It is useful 
to note that a similar error term would exist in Equation 8 for a retrieval algorithm that 
did not correct for w[ in Equation 7. Measurement errors in Vj and V 2 would generate 
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additional error terms. (u 2 + u 2 )r*/ be a ff ec t £ d by the error term because turbulent 
eddies cause the vertical velocity to vary rapidly in time and space, resulting in w 2 * w, . 
As a result, (U 2 + u 2 ) Rer - can be expected to be more variable than (u 2 +u 2 ) True . 
However, (U 2 +u 2 ) Rel will be unbiased with respect to (U 2 + u 2 ) True if w 2 - w, is zero 
and if variations in the error term are uncorrelated with perturbations in (U 2 ) True : 

On the other hand, (u° \ el may be systematically biased if positive peaks in the error 
term coincide with peak or near-peak values in (u 2 + u 2 ) 7n ,„ • The probability of such 
coincidence would increase under one or more of the following three conditions: 1) As 

the averaging interval becomes long, 2) As the time-scale of w variations becomes short 

compared to the time-interval between observations, or 3) As the distance between the 
beams becomes large compared to the spatial scale of the turbulent eddies. Some 
quantitative insights into these potential errors of (uf ) Re( can be obtained by use of a 
statistical model, as described below, to simulate Doppler profiler observations and the 
resulting (u 2 +u 2 ) Rel . 

3. Statistical Modeling of Peak Horizontal Velocity Bias 

The idealized profiler concept introduced in Section 2 will be used here to obtain 

quantitative insights into the statistics of and i U n L b y employing analytical 

properties of the Extreme Value Type-I distribution for maxima (the Gumbel 
distribution). 
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A Gumbel distribution can be obtained by generating random samples of size n from 
a normal distribution, then extracting the maximum value from each sample and 
repeating the process ad infinitum (Coles 2001). The maxima will have a Gumbel 

distribution. In the present case, (t/ 2 + u 2 ) True is assumed to be normally distributed 

about TT 2 , with a standard deviation of ct (Mitsuta and Tsukamoto 1989). In a thought 
experiment we imagine that from each random sample of size n from the population of 
{u 2 + u 2 \ me , the peak, or maximum, value is selected and used to create a population of 
{u®\ ne . The distribution of (u* \ ne will be Gumbel in form with scale and location 
parameters X,x r ue and ^rrue, respectively. At the same time a sample of {u 2 + u 2 ) Re( is 
generated consisting of {lJ 2 + u 2 ) Tnie plus the random error term as described by Equation 
8. The error term comes from normal distributions of wjandw, with zero means and 
equal variances. A corresponding distribution of {u n ) Re( is generated by taking the 
maximum (f^ + « 2 ) Re( for each sample of size n. When the error term is normally 
distributed and independent of (u 2 + u 2 ) True , the resulting distribution of will 

also be Gumbel in form with parameters X. Re t and £ Re t- The equations below describe 
how X True , £ True , X Ket and £ Ret depend on n, U 2 , cj, and parameters of the error term. 
The error term will have a variance that depends on properties of (w 2 -w,) and the 

amplifying cot(©) factor. 

The Gumbel PDF is given by 

G( £7® ) = 1/ X • exp{-(£/®-^)/ X } • exp[-exp{-(£/® -£)/ X }].(9) 
The average value and variance of the maxima, £7® in this case, are given by 
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(la®) = ^ + X • 0.5772 and 
Var [U®]=X 2 -n 2 /6, 


( 10 ) 


( 11 ) 

respectively. The value 0.5772 is an approximation of Euler’s constant (Abramowitz and 
Stegun 1965), hereafter Eu. 

When U 2 and CTtrue are the mean and standard deviation of (u 2 + u 2 ) Tnie , the 
parameters of the resulting Gumbel PDF of (u® \ ne are 


^ True ^2 ^True 

and 

(12) 

= b n * CJxrue? 


(13) 


where, a„ and b„; known as normalizing constants, are dependent on the sample size and 
are given by (T. Rolf Turner, University of New Brunswick, Canada, personal 
communication): 

a„ = sqrt[2ln(n)] - [ln(4rc) + ln(ln(n))]/[2 sqrt(21n(n))] (14) 

b n = l/sqrt[21n(n)]. (15) 

The average value of {u®\ ne can then be written in terms of parameters of the 
underlying normal distribution and the normalizing constants by 

(LA® J^~ e = U 2 + O' True ‘ (Sn'I’Eubn). (16) 

The (a n +E u b n ) factor, hereafter G c , is weakly dependent on n, changing from 2.56 to 
2.92 as n goes from 100 to 300. For example, consider U 2 = 10, CT T rue = 2.0, and n=100. 

We find = 15.1. 
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Gumbel distribution parameters for (t/® ) ReI are made by making use of Equation 8 in 

Section 2b. The variance of ^T 2 + u 2 ) Re , will be affected by the error term, due to 

variance of (w 2 -wj) and the covariance of ( U 2 +u 2 ) with (w 2 -w ] ). The problem of 

solving for the variance of (u 2 + u 2 ) Re , and for the Gumbel parameters of (u®\ et can be 

simplified by assuming that the covariance between ( U 2 + u 2 ) and ( w 2 - w i ) is zero. That 

is, the difference in vertical velocities over the bi- and b 2 -beams is assumed to be 
uncorrelated with perturbations in the horizontal velocity. This assumption will be made, 

causing the variance of {u 2 +w 2 ) Ref to be greater than that of (f/ 2 + w 2 ) rrue . It is the 
increased variance of (t/ 2 + u 2 that results in a positive bias in (f/f • 

The additional variance in ( U 2 + u 2 ) Rel will be dependent on the variance of the 

difference in vertical velocities over. the two beams, (w 2 -w,). The variance of this 
difference will be referred to as the variance of the residual, Var[Resid] because it could 
be zero if the vertical velocities over the two beams were identical. In the case where w 2 
is known the retrieval equation will not include the error term. 

The variance of + u 2 ) Re , , Var[U Ret ], is then obtained by summing the variance of 

frj~ 2 + u 2 ) Tme , or (jj me , with the product of Var[Resid] and the square of the amplifying 

cotangent term: 

Var[U R et] = cT^ ne + Var[Resid]- Cot 2 (0). (17) 

For the case where Var[Resid] is (1.0)“, 0 = 15°, and the variance of ( U 2 +w 2 )rrue 
(2.0) 2 , Var[U Ret ], becomes 17.9, a substantial increase. 
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The average value for the (t/® ) Re( can now be calculated from 


WL = U 1 ( 18 > 

where <7 Ret = {Var[U Ret j} 1/2 . From Equation 18 we find (t/f = 20.8, using U 2 =10, 

a 2 = 4.0, and n=100, as above, is -38% higher than (U ® ) rn/e , which was found 

to be 15.1. 

Figure 2 shows modeled probability density functions (pdfs) of {lJ 2 + u 2 ) Tnie , 
and (f/®L for the case outlined above. The {u®l ne and (u® distributions 

are shifted significantly to the right of the underlying distribution of (u 2 + u 2 ) True , 
because the theoretical peak value distributions were determined for sample sizes of 100. 
It is also evident that the distribution of (u® \ 6l is somewhat wider than that of (u® . 

For the case where 0 = 30°, Var[U Ret ] would be reduced to 7.0 and the bias would be 
reduced to 11%. The bias decreases to zero as 0 approaches 90°, as expected. 

The influence of the variance of ( w 2 - w, ) in enhancing (u® ) Re( above (u® ) rrue can be 
visualized by graphing the percent bias of (f/® j^ e< in a two-dimensional, normalized 
coordinate system of a Tru JU 2 versus CT Res id IU 2 . Figure 3 shows contours of the 
percent-bias, 

%Bias = 100*(^L - KL VWL. ( A12 ) 

for the case where n=100, and theta = 15 degrees. The greatest bias is found when the 
residual variance in vertical velocity is large compared to the variance in horizontal 
velocity. The bias becomes small when the residual variance in vertical velocity is small. 
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Figure 4 shows contours of the percent bias for the case where n=100 and 0 = 30°. 
For fixed values of a Resid /U 2 and c r /U 2 , the percent bias with 0 = 30° is less than half 
that for 0 = 15° (Fig. 3) because of the strong © dependence of the cot 2 (0) factor. 

4. Summary and Conclusions 

Numerous previous studies have shown that Doppler profilers are capable of 
providing accurate retrievals of average horizontal velocities in the atmosphere and 
underwater. However, estimates of peak horizontal velocities from collections of 
instantaneous retrievals are susceptible to a positive bias, due to turbulent vertical 
motions in the medium of interest. 

An idealized Doppler profiler configuration was combined with a statistical model 
of turbulent motions and a simple retrieval algorithm to illustrate the nature of the bias 
and its magnitude. The results revealed that unresolved vertical motions contaminate the 
instantaneous retrievals because of limitations inherent in the beam configuration of 
typical profilers. The retrieved instantaneous horizontal velocities are more variable than 
the true instantaneous velocities, resulting in a positive bias when the peak retrieved 
value is chosen. 

The bias in peak horizontal velocity can be characterized in terms of turbulent 
properties of the flow. Errors in the measurement of the vertical velocity would also 
contribute to the bias. These findings suggest that an average correction could be applied 
to the retrieved peak values. However, the correction would be statistical in nature and 
would not necessarily improve the precision of individual retrieved peak values. 


13 



APPENDIX A 


List of Symbols 


a n 

Gumbel normalization constant 

bn, 

Gumbel normalization constant 

bj-beam 

Vertical beam 

b2-beam 

Oblique beam 

Di 

Distance along bi-beam axis 

d 2 

Distance along b 2 -beam axis 

Eu 

Euler’s constant 

G c 

Gumbel constant (a n +b n -Eu) 

G(x) 

Gumbel probability density function 

H 

Height above the profiler 

n 

Sample size 

P 

Location of profiler 

Resid 

Residual quantity 

R 

Radial velocity 

r 

Perturbation R 

U 

Horizontal velocity 

u 

Perturbation horizontal velocity 

Var[x] - 

Variance of the variable x 

V, 

Radial velocity along bi-beam axis 

V 2 

Radial velocity along b 2 -beam axis 

w 

Vertical velocity 
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w 


0 


5 


1 

2 


Ret 


True 


Perturbation vertical velocity 

-j 

Angle between vertical and oblique beams 
Scale parameter for Gumbel distribution 

Location parameter for Gumbel distribution 
Standard deviation 

Overbar indicates ensemble average value 
Prime indicates perturbation value 
Indicates a peak value 

Subscript 1 indicates vertical beam 
Subscript 2 indicates oblique beam 
Subscript Ret indicates a retrieved value 
Subscript True indicates a true value 
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Figure Captions 


Figure 1 . Schematic of an idealized 2-axis Doppler profiler system 
Figure 2. Modeled pdfs of horizontal velocities illustrating Gumbel distributions of true 
and retrieved peaks derived from a background distribution that is Gaussian. 
Distribution parameters have been chosen to give an average background velocity of 
10, an average value of 15.1 for the true peaks and a bias of 38% for the retrieved 
peaks. The standard deviation of true horizontal velocities is 2.0 and the standard 
deviation of the residual vertical velocity is 1 .0. 

Fig. 3. The percent bias in \Uf j^ e( as a function of crime / U and CJResid !U , for n=100 

and 0 = 15°, based on the analytical models for horizontal and vertical velocity 
presented in Section 3. 

Fig. 4. As in Fig.3, but for © = 30°. 
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p 

Figure 1. Schematic of an idealized 2-axis Doppler profiler system. 


21 



Modeled PDFs of Horizontal Velocities: 
Background, True Peaks and Retrieved Peaks 
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Figure 2. Modeled pdfs of horizontal velocities illustrating Gumbel distributions 
of true and retrieved peaks derived from a background distribution that is 
Gaussian. Distribution parameters have been chosen to give an average 
background velocity of 10, an average value of 15.1 for the true peaks and a bias 
of 38% for the retrieved peaks. The standard deviation of true horizontal 
velocities is 2.0 and the standard deviation of the residual vertical velocity is 1.0. 
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% Bias: Average Retrieved (U*) : n=100, theta= 15 



Figure 3. The percent bias in (£/®) Re , as a function of csj ru JU and aR es id /U , for 

n=100 and © = 15°, based on the analytical models for horizontal and vertical velocity 
presented in Section 3. 
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